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Abstract. Dimuon production has been studied in a series of fixed-target experiments at Fermilab 
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during the last two decades. Highlights from these experiments, together wih recent results from the 
Fermilab E866 experiment, are presented. Future prospects for studying the parton distributions in 
the nucleons and nuclei using dimuon production at Fermilab and J-PARC are also discussed. 
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INTRODUCTION 

X\ 
i ■ While our current knowledge on the partonic substructures of nucleons and nuclei is 

q ■ obtained mostly from deep-inelastic scattering experiments, the dimuon production ex- 

periments also play an important and unique role in providing complementary informa- 
tion. During the last two decades, a series of fixed-target dimuon production experiments 
(E772, E789, E866) have been carried out using 800 GeV/c proton beam at Fermilab. 
At 800 GeV/c, the dimuon data contain Drell-Yan continuum up to dimuon mass of 
~ 15 GeV as well as quarkonium productions (SP¥, *¥', and T resonances). The Drell- 
Yan process and quarkonium productions often provide complementary information, 
since Drell-Yan is an electromagnetic process via quark-antiquark annihilation while 
the quarkonium production is a strong interaction process dominated by gluon-gluon 
fusion at this beam energy. 

The Fermilab dimuon experiments covers a broad range of physics topics. The Drell- 
Yan data have provided informations on the antiquark distributions in the nucleons yi 
l2Ll3i|4{] and nuclei [5, 6]. These results showed the surprising results that the antiquark 
distributions in the nuclei are not enhanced [5, 6], contrary to the predictions of models 
which explain the EMC effect in term of nuclear enhancement of exchanged mesons. 
Moreover, the Drell-Yan cross section ratios p + d/p + p clearly establish the flavor 
asymmetry of the d and u distributions in the proton, and they map out the x-dependence 
of this asymmetry D2l,[3|,|4J]. Pronounced nuclear dependences of quarkonium productions 
have been observed for J/ 1 !', *¥', and T resonances Q7|» ls|, |9L llO(|. It was found that these 
nuclear effects scale with the kinematic variable xf rather than X2 (the Bjorken-x of the 
parton in the nulceus), suggesting the importance of initial- and final-state interactions. 
A striking behavior of the nuclear dependence as a function of pr, reminiscent of the 
Cronin effect, was also observed. The nuclear Drell-Yan cross sections also exhibit 
xf as well as pr dependences [a luLllJ, LLSQ, which are weaker than the quarkonium 
nuclear dependences but can provide information on the energy loss of quarks traversing 
the nucleus [14]. The differential cross sections for Drell-Yan [15], charmonium UM, 
and bottomonium Ill7n productions have also been reported. In addition, the decay 



angular distributions for Drell-Yan till Il8ll . JW UM, and T resonances [20] have 



been measured. Information on D and B meson productions and decays has also been 



extracted through either the dimuon measurement 112 ll . |22| . 12311 or the open-aperture 



dihadron measurement D24J]. Several review articles covering some of these results are 
available rillLl25l 12611 . In this article, we will focus on the recent results from experiment 



E866 and future prospect of dimuon experiments at Fermilab and J-PARC. 

RECENT RESULTS FROM E866 

Angular distributions of proton-induced Drell-Yan 

Despite the success of perturbative QCD in describing the Drell-Yan cross sections, 
it remains a challenge to understand the angular distributions of the Drell-Yan process. 
Assuming dominance of the single-photon process, a general expression for the Drell- 
Yan angular distribution is I127n 

do . 7 v 7 

— oc 1 +Acos + usin20cosd> + — sin 0cos20, (1) 

ail 2 

where and denote the polar and azimuthal angle, respectively, of the Z + in the 
dilepton rest frame. In the "naive" Drell-Yan model, where the transverse momentum 
of the quark is ignored and no gluon emission is considered, A = 1 and \x — V = 
are obtained. QCD effects I128il and non-zero intrinsic transverse momentum of the 
quarks [|29|l can both lead to X 7^ 1 and /1, v 7^ 0. However, X and V should still satisfy the 
relation 1 — X = 2v [27]. This so-called Lam-Tung relation, obtained as a consequence 
of the spin- 1/2 nature of the quarks, is analogous to the Callan-Gross relation [30] in 
deep-inelastic scattering. 

The first measurement of the Drell-Yan angular distribution was performed by the 



NA10 Collaboration for it' + W with the highest statistics at 194 GeV/c |31|,|32Q. The 
cos 20 angular dependences showed a sizable v, increasing with dimuon transverse 
momentum (pj) and reaching a value of?«0.3atpr = 2.5 GeV/c. The Fermilab E615 
Collaboration subsequently performed a measurement of %~ + W Drell-Yan production 



at 252 GeV/c with broad coverage in the decay angle (330. The E615 data showed that 
the Lam-Tung relation, 2v = 1 — X, is clearly violated. 

The NA10 and E615 results on the Drell-Yan angular distributions strongly sug- 
gest that new effects beyond conventional perturbative QCD are present. Brandenburg, 
Nachtmann and Mirke suggested that a factorization-breaking QCD vacuum may lead 
to a correlation between the transverse spin of the antiquark in the pion and that of the 
quark in the nucleon [34]. This would result in a non-zero cos 20 angular dependence 
consistent with the data. Several authors have also considered higher-twist effects from 



quark-antiquark binding in pions |35|,|36J]. However, the model is strictly applicable only 
in the x n — ► 1 region, while the NA10 and E615 data exhibit nonperturbative effects over 
a much broader kinematic region. 

More recently, Boer pointed out [37] that the cos 20 angular dependences observed in 
NA10 and E615 could be due to the /^-dependent parton distribution function h\. This 
so-called Boer-Mulders function [38] is an example of a novel type of fc^-dependent 
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FIGURE 1. Parameters A,/x, V and2v— (1 — A) vs. pj in the Collins-Soper frame. Solid circles are for 
E866 p + d at 800 GeV/c, crosses are for NA10 n~ + W at 194 GeV/c, and diamonds are E615 n~ + W 
at 252 GeV/c. The error bars include the statistical uncertainties only. 



parton distribution function, and it characterizes the correlation of a quark's transverse 
spin and its transverse momentum, kj, in an unpolarized nucleon. 

To shed additional light on the origins of the NA10 and E615 Drell-Yan angular 
distributions, we have analyzed p + p and p + d Drell-Yan angular distribution data 
at 800 GeV/c from Fermilab E866. There has been no report on the azimuthal angular 
distributions for proton-induced Drell-Yan - all measurements so far have been for polar 
angular distributions Illll.ll9l. l20n. Moreover, proton-induced Drell-Yan data provide a 
test of theoretical models. For example, the cos 20 dependence is expected to be much 
reduced in proton-induced Drell-Yan if the underlying mechanism involves the Boer- 
Mulders functions. This is due to the expectation that the Boer-Mulders functions are 
small for the sea-quarks. However, if the QCD vacuum effect 113411 is the origin of the 
cos 20 angular dependence, then the azimuthal behavior of proton-induced Drell-Yan 
should be similar to that of pion-induced Drell-Yan. Finally, the validity of the Lam- 
Tung relation has never been tested for proton-induced Drell-Yan, and the present study 
provides a first test. 
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FIGURE 2. Parameter v vs. pj in the Collins-Soper frame for several Drell-Yan measurements. Fits to 
the data using Eq. 3 and Mq = 2.4 GeV/c 2 are also shown. 



The Fermilab E866 experiment was performed using the upgraded Meson-East mag- 
netic pair spectrometer. Details of the experimental setup have been described else- 
where [0, SH]. An 800 GeV/c primary proton beam with up to 2 x 10 12 protons per 
20-second beam spill was incident upon 50.8 cm long cylindrical stainless steel target 
flask containing liquid deuterium. The detector system consisted of four tracking stations 
and a momentum analyzing magnet. From the momenta of the /i + and [i~, kinematic 
variables of the dimuons (xf .m^ : pt) were readily reconstructed. The muon angles 6 
and in the Collins-Soper frame [39] were also calculated. To remove the quarkonium 



background, only events with 4.5 < mnu, < 9 GeV/c or man > 10.7 GeV/c were an- 
alyzed. A total of ~54,000 p + p and ~1 18,000 p + d Drell-Yan events covering the 
decay angular range —0.5 < cos < 0.5 and — % < (j> < % remain. 

Figure 1 shows the results [18] on the angular distribution parameters A,/i, and V 
vs. pj. To extract these parameters, the Drell-Yan data were grouped into 5 bins in 
cos 6 and 8 bins in for each pj bin. A least-squares fit to the data using Eq. 1 to 
describe the angular distribution was performed. Only statistical errors are shown in 
Fig. 1. For comparison with the p + d Drell-Yan data, the NA10 K~ + W data at 194 
GeV/c and the E615 %~ + W data at 252 GeV/c are also shown in Fig. 1. To test the 
validity of the Lam- Tung relation, also shown in Fig. 1 is the quantity, 2 v — ( 1 — X ) , for 



all three experiments. For p + d at 800 GeV/c, Fig. 1 shows that A is consistent with 
1, in agreement with previous studies |[LiLlL9|,|2Q1], while /i and V deviate only slightly 
from zero. This is in contrast to the pion-induced Drell-Yan results, in which much larger 
values of V are found. It is also interesting to note that while E615 clearly establishes the 
violation of the Lam- Tung relation, the NA10 and the p + d data are largely consistent 
with the Lam-Tung relation. 

The p + d results put constraints on theoretical models that predict large cos 20 
dependence originating from QCD vacuum effects. They also suggest that the Boer- 
Mulders function h\ for sea quarks is significantly smaller than for valence quarks. A 
recent analysis [40] of the p + d cos 2(j) data showed that the sea-quark Boer-Mulders 
functions are indeed smaller by a factor ~ 5 than the valence-quark Boer-Mulders 
functions. This analysis also indicated that the E866 p + d data are consistent with the 
u and d Boer-Mulders functions having the same signs. However, the p + d data alone 
can not provide an unambiguous determination of the flavor dependence of the Boer- 
Mulders functions. A comparison of the p + p and p + d data would further constrain 
the flavor dependence of the Boer-Mulders functions. 

Figure 2 shows the preliminary result of V for p + p reaction at 800 GeV/c. Also 
shown in Fig. 2 are the E866 p + d and the NA10 n~ + W data at 194 GeV/c. The data 
are fitted with the expression suggested by Boer [37]: 

v = 16ki , C 9N0 , (2) 

where K"i is proportional to the product of the hj- functions for the projectile and the 
target, and Mq is a constant fitting parameter. Boer obtained K"i = 0.47 ±0.14 and 
M c = 2.4 ± 0.5 GeV/c 2 for fitting the NA10 data, as shown in Fig. 2. A fit to the E866 
p + p and p + d data for M c = 2.4 GeV/c 2 yields K\ = 0.21 ±0.055 and K X = 0. 1 1 ±0.04, 
respectively. These data should provide additional information on the origins of the 
azimuthal angular dependence of Drell-Yan, as well as the flavor dependence of the 
Boer- Mulders functions. 



T production for p + p and p + d interactions 

In the CERN NA51 lEH] and Fermilab E866 10, fill] experiments on proton-induced 
dimuon production, a striking difference was observed for the Drell-Yan cross sec- 
tions between p+p and p + d. As the underlying mechanism for the Drell-Yan pro- 
cess involves quark-antiquark annihilation, this difference has been attributed to the 
asymmetry between the up and down sea quark distributions in the proton. From the 
o(p + d)Dy/2o(p + p)dy ratios the Biorken-jc dependence of the sea-quark d/u flavor 
asymmetry has been extracted p.l3l.H. l4l|l . 

The Fermilab E866 dimuon experiment also recorded a large number of T — > ll + fi 
events. Unlike the electromagnetic Drell-Yan process, quarkonium production is a strong 



interaction dominated by the subprocess of gluon-gluon fusion at this beam energy [|42 . 
14311 . Therefore, the quarkonium production cross sections are primarily sensitive to 
the gluon distributions in the colliding hadrons. The T production ratio, o(p + d — * 
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FIGURE 3. The E866 o(p + d) /2a (p + p) cross section ratios for T resonances as a function of x%. 
The corresponding ratios for Drell-Yan cross sections are also shown. The error bars are statistical only. 



T)/2a(p + p — > T), is expected to probe the gluon content in the neutron relative to that 
in the proton [44] . While it is generally assumed that the gluon distributions in the proton 
and neutron are identical, this assumption is not based on any fundamental symmetry and 
has not been tested experimentally. A possible mechanism for generating different gluon 
distributions in the proton and neutron, as pointed out by Piller and Thomas [44], is the 
violation of charge symmetry in the quark and antiquark distributions in the nucleons. 
A precise measurement of the o(p + d —> T)/2o(p + p —> Y) ratios would provide a 
constraint on the asymmetry of gluon distribution in the proton versus that in the neutron. 
The o(p + d)/2o(p + p) ratios for T( 15 + 25 + 35) production are shown in Fig. 3 as 
a function of x%. Most of the systematic errors cancel for these ratios, with a remaining 
«1% error from the rate dependence and target compositions [4]. Figure 3 shows that 
these ratios are consistent with unity, in striking contrast to the corresponding values D4(] 
for the Drell-Yan process, also shown in Fig. 3. The difference between the Drell- 
Yan and the T cross section ratios clearly reflect the different underlying mechanisms 
in these two processes. The Drell-Yan process, dominated by the u — u annihilation 
subprocess, leads to the relation <j(p + d)Dy/2a(p + p)dy ~ §(1 + u n (x2) /tipfa)) = 
i(l + dp(x2)/up(x-i)), where q p ^ refers to the q distribution in the proton and neutron, 
respectively. For T production, the dominance of the gluon-gluon fusion subprocess at 
this beam energy implies that o(p + d — ► T)/2o(p + p — ► Y) « ^(1 + g n {x2) / 'gp(x2)) ■ 
Figure 3 shows that the gluon distributions in the proton (g p ) and neutron (g n ) are very 
similar over the X2 range 0.09 < X2 < 0.25. The overall o(p + d — > T)/2o(p + p — > Y) 
ratio, integrated over the measured kinematic range, is 0.984+0. 026(stat.) ±0.01(syst). 
The Y data indicate that the gluon distributions in the proton and neutron are very 
similar. These results are consistent with no charge symmetry breaking effect in the 
gluon distributions. 



FUTURE PROSPECTS AT FRMILAB AND J-PARC 

Future fixed-target dimuon experiments have been proposed at the 120 GeV Fermilab 
Main Injector (FMI) and the 50 GeV J-PARC facilities. The Fermilab proposal Q, 
E906, has been approved and is expected to start data-taking around 201 1. Two dimuon 
proposals (P04 [46] and P24 [47]) have also been submitted to the J-PARC for approval. 
The lower beam energies at FMI and J-PARC present opportunities for extending the 
d/u and the nuclear antiquark distribution measurements to larger x (x > 0.25). For 
given values of x\ and x%, the Drell-Yan cross section is proportional to \/s, hence 
a gain of ~ 16 times in the Drell-Yan cross sections can be obtained at the J-PARC 
energy of 50 GeV. Since the perturbative process gives a symmetric d/u while non- 
perturbative processes are necessary to generate an asymmetric d/u sea, it would be 
very important to extend the Drell-Yan measurements to kinematic regimes beyond the 
current limits. Another advantage of lower beam energies is that a much more sensitive 
study of the partonic energy loss in nuclei could be carried out using the Drell-Yan 



nuclear dependence [14]. 

The dimuon physics program at J-PARC is proposed to be carried out in several 
stages. Since 30 GeV proton beam will be available at the initial phase of J-PARC, the 
first measurements will focus on J /*¥ production at 30 GeV. This will be followed by 
measurements of Drell-Yan and quarkonium production at 50 GeV after the beam energy 
is upgraded to 50 GeV. Experiments using polarized target could already be performed 
with unpolarized beams. When polarized proton beam becomes available at J-PARC, a 
rich and unique program on spin physics could also be pursued at J-PARC using the 
dimuon spectrometer. 

An important feature of J '/*¥ production using 30 or 50 GeV proton beam is the 
dominance of the quark-antiquark annihilation subprocess. This is in striking contrast 
to Jp9 production at 800 GeV (Fermilan E866) or at 120 GeV (Fermilab E906), where 
the gluon-gluon fusion is the dominant process. This suggests an exciting opportunity to 
use J '/*P production at J-PARC as an alternative method to probe antiquark distribution. 

With the possibility to accelerate polarized proton beams at J-PARC, the spin struc- 
ture of the proton can also be investigated with the proposed dimuon experiments. In 
particular, polarzied Drell-Yan process with polarized beam and/or polarized target at 
J-PARC would allow a unique program on spin physics complementary to polarized 
DIS experiments and the RHIC-Spin program. Specific physics topics include the mea- 
surements of T-odd Boer-Mulders distribution function in unpolarized Drell-Yan, the 
extraction of T-odd Sivers distribution functions in singly transversely polarized Drell- 
Yan, the helicity distribution of antiqaurks in doubly longitudinally polarized Drell-Yan, 
and the transversity distribution in doubly transversely polarized Drell-Yan. It is worth 
noting that polarized Drell-Yan is one of the major physics program at the GSI Polar- 
ized Antiproton Experiment (PAX). The RHIC-Spin program will likely provide the first 
results on polarized Drell-Yan. However, the high luminosity and the broad kinematic 
coverage for the large-x region at J-PARC would allow some unique measurements to 
be performed in the J-PARC dimuon experiments. 

I am very grateful to my many collaborators on the E772, E789, and E866 experiments 
at Fermilab. I would also acknowledge the collaboration with Dr. Shin'ya Sawada and 
Dr. Yuji Goto on the J-PARC dimuon proposals. 
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